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Abstract. State of the art research and treatment of biological tissues require accurate and 
efficient methods for describing their mechanical properties. Indeed, micromechanics moti- 
vated approaches provide a systematic method for elevating relevant data from the microscopic 
level to the macroscopic one. In this work the mechanical responses of hyperelastic tissues 
with one and two families of collagen fibers are analyzed by application of a new variational 
estimate accounting for their histology and the behaviors of their constituents. The resulting, 
close form expressions, are used to determine the overall response of the wall of a healthy hu- 
man coronary artery. To demonstrate the accuracy of the proposed method these predictions are 
compared with corresponding 3-D finite element simulations of a periodic unit cell of the tissue 
with two families of fibers. Throughout, the analytical predictions for the highly nonlinear and 
anisotropic tissue are in agreement with the numerical simulations. 
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1. Introduction 

Arteries and other blood vessels, tendons, ligaments, skin, cornea, cartilage, intervertebral 
discs and the gut, among others, are deformable collagenous tissues whose physiological func- 
tioning crucially depends on their mechanical properties. Moreover, it is well known that 
the mechanical environment experienced by these tissues play a critical role in the origin and 
progression of diseases. In particular, with respect to the cardiovascular system: stroke, aortic 
aneurisms and dissections, acute myocardial infraction (MI) and congestive heart failure (CHF) 
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as well as restenosis following intravascular stenting and other cardiovascular diseases and in- 
juries, result in significant morbidity and mortality worldwide. Accordingly, many methods 
of clinical treatment and surgical interventions largely rely on mechanical effects (Humphrey, 
2002). Over the last decade, the unique mechanical behavior of soft biological tissues and 
particularly that of the arterial wall has been receiving considerable attention in the literature. 
Despite the significant progress, there remains a pressing need for a better understanding of 
the fundamental role of mechanics in vascular physiology and pathophysiology, that will even- 
tually enable the prediction of the mechano- stimulated vascular adaptation (Humphrey, 2003; 
Holzapfel and Ogden, 2010). 

On top of the complications related to the biological aspects, even the characterization of the 
purely mechanical factors governing the response of soft connective tissues is not a straight- 
forward procedure and requires a reliable constitutive model. The main constituents of soft 
tissues are: fibroblasts, the extra-fibrillar matrix (EFM) which includes the elastin, and the col- 
lagen fibers. From a mechanical viewpoint, the most important load-bearing component is the 
collagen fibers network. The function, strength, integrity and stability of soft tissues are main- 
tained by the structural arrangement of the collagen fibers. Thus, changes in collagen content, 
orientation, fiber type or fiber thickness, all have a major impact on the overall responses of the 
tissues. Many of the aspects related to the collagen fibers are discussed in Fratzl (2008). 

The histology of collagenous tissues evolves in a way that optimizes their overall mechanical 
functionality. Specifically, the mechanical behavior of the artery depends on its anatomical site. 
For instance, Taber and Humphrey (2001) note that the properties of the vascular tissues vary 
across the arterial wall. Also, it has been demonstrated that the distribution of the angle of the 
collagen fibers network vary through the thickness of the artery, particularly at bifurcation sites 
(Holzapfel et al., 2002; Hariton et al., 2007b). 

Today it is widely accepted that changes in the ambient mechanical loads acting on the tissue 
are followed by a remarkable continuous adaptation of the tissue. The versatility of collagen 
as a building material is mainly thanks to modification of its structure. Collagen adaptation 
to loading both in healthy and diseased tissues involves two main processes: remodeling and 
growth. Particularly, the fibers may grow and reorient in order to meet new mechanical de- 
mands. Growth is defined as the change in concentration due to mass transport. Remodeling is 
considered as reorganization of existing fibers at constant mass (Ambrosi et al., 2011). 

The list of examples of vascular growth and remodeling in normal development, ageing, 
disease progression, functional adaptations to altered environments and responses to injury or 
clinical treatment is endless (e.g., Ambrosi et al., 2011). Hypertension is an example of a 
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disease that causes and is caused by alternations in the mechanical environment of the arterial 
wall. In hypertension, the sustained increase in blood pressure may result in a thickening 
of the wall as a consequence of collagen deposition and associated restoration of the hoop 
stress to its normal value. There are also changes in blood flow-induced wall shear stress that 
effect the artery wall such as enlargement of lumen diameter due to increased wall shear stress. 
Obviously, during the development of hypertension, marked changes in the microstructure and 
functionality of blood vessels take place. Growth and remodeling processes of collagen fibers 
serve to restore the stresses to their homeostatic levels. These processes attempt to recover and 
maintain the preferred form and function and are fundamental for adaptation to altered blood 
flow (Humphrey, 2008; Eberth et al., 2010). In the context of tissue engineering, Nerem and 
Seliktar (2001) emphasized the strong influence of mechanical loading on the remodeling of 
the collagen network. 

It is evident that in order to accurately capture the stress-strain response of the heteroge- 
neous, anisotropic, and highly nonlinear tissue with its characteristic stiffening at large strains, 
it is necessary to develop a tool that enables to account for the specific contributions of the 
significant constituents to its overall mechanics. Moreover, to conduct remodeling analyses 
that provide further insight into the mechanical behavior of the phases composing the tissue, 
structurally-based constitutive models that incorporate collagen fiber alignments, are needed. 
The motivation for our research is best summarized in Humphrey (2009) who states that one 
of the conspicuous shortcomings that have remained is "the need to account explicitly for sep- 
arate contributions by the diverse structurally significant constituents within the arterial wall". 
Nonetheless, while phenomenological models that treat the arterial walls and other soft col- 
lagenous tissues as homogeneous materials were developed in recent years (e.g., Fung, 1967; 
Lanir, 1983; Humphrey and Yin, 1987; Holzapfel et al., 2000, among others), to date only a 
few micromechanics-motivated studies exist in the literature (e.g., deBotton and Shmuel, 2009; 
Chen et al.,2011). 

The aim of this work is to introduce a novel comprehensive micromechanical-based model 
for the determination of the macroscopic response of soft tissues with two families of fibers in 
terms of the properties of the constituting phases, their spatial distribution and the interactions 
between them in a consistent way. At the present stage of the study we restrict our attention 
to passive response of tissues which are not pre-stressed in their load free configuration, and 
neglect aspects such as fiber-fiber and fiber-matrix cross-links (e.g., Avery and Bailey, 2008). 
Nonetheless, adaptation of this approach allows to reflect features of the arteries geometry 
and morphology, in particular the orientations distinguished by the collagen fiber structure and 
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variations in their content. Among other aspects this approach will also allow to account for 
changes in the composition of the constituents, thickening of the fibers, alternation in their 
density, deterioration of the EFM due to disease (e.g., de Figueiredo Borges et al., 2008) or 
variations in nutrition habits. Indeed, the advantages of being able to extract the overall be- 
havior of the tissue from the individual behaviors of the phases and their histology are clear. 
In particular, if fair estimates for the fiber and the extra-fibrillar matrix behaviors are known, 
this technique will play a key role in providing exact and realistic predictions for the mechan- 
ical properties of both healthy and diseased tissues that can be deduced from patient- specific 
non-invasive inspection techniques. 

Consistent micromechanics analyses are commonly based on homogenization theories. While 
numerous works adopting this approach exist in the limit of infinitesimal deformation elastic- 
ity, for finitely deforming materials only a few works can be found in the literature. The 
foundations of the homogenization theory for geometrically nonlinear materials were outlined 
40 years ago by Hill (1972) and Ogden (1974), however, analytical results for fibrous materials 
that are based on this theory become available only in the last decade. Rigorous predictions 
and estimates for hyperelastic materials with one family of fibers were introduced by deBotton 
and Hariton (2006); deBotton et al. (2006); Shmuel and deBotton (2010) and Lopez-Pamies 
and Idiart (2010), and corresponding results for materials with two families of fibers were de- 
termined by deBotton and Shmuel (2009). Estimates for fiber composites that are based on the 
linear comparison method of Ponte Castaneda and Tiberio (2000) were determined by Brun 
et al. (2007) and Agoras et al. (2009). Chen et al. (2011) followed this method to estimate 
the behavior of materials with two fiber families. The results presented in this work are based 
on the nonlinear comparison (NLC) variational estimate introduced recently by deBotton and 
Shmuel (2010). This method allows to extend available estimates for the macroscopic behav- 
ior of hyperleastic heterogeneous materials to families of hyperelastic materials with different 
constitutive behaviors of the phases. 

To demonstrate the validity of the estimates obtained by application of the NLC variational 
method we compare its predictions with corresponding finite element (FE) simulations. In this 
regard we note that due to the expected advantage of the micromechanics-motivated approach, 
the microsphere formulation, which is a finite element based technique that accounts for the 
local microstructure of the material, was recently applied to model the response and remodeling 
in soft collagenous tissues (e.g., Menzel and Waffenschmidt, 2009; Alastrue et al., 2009). Here, 
we numerically tackle the homogenization problem by considering a 3-D periodic specimen 
whose representative unit cell contains two families of fibers. We note that this approach 
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is inefficient for large models, but in this work it allows to compare the predictions of the 
analytical estimates under uniform loading conditions. 

The work is structured as follows: In the following section we briefly revisit the required 
background including homogenization theory in finite deformation elasticity. We also recall 
the variational statement that enables to determine the analytical estimates for the behavior 
of the collagenous tissue as well as the complementary finite element procedure for periodic 
media. In section 3 we apply the two techniques to a specific tissue, namely a healthy hu- 
man coronary artery. This is accomplished by appropriate choices of the behaviors of the 
constituents composing the artery and its histology. Finally, predictions for the response of 
the tissue to different loadings according to the analytical and the numerical formulations are 
determined and compared. A short summary concludes the work. 

2. THEORETICAL BACKGROUND 

The Cartesian position vector of a material point in a reference configuration of a body 
is X, and its position vector in the deformed configuration SB is x. The deformation of the body 
is characterized by the mapping 

x = *(X). (1) 

The deformation gradient is 

We assume that the deformation is invertible and hence F is non-singular, accordingly 

/ = detF ^ 0. (3) 

Physically, / is the volume ratio between the volumes of an element in the deformed and the 
reference configurations, and hence J = ^ >0. 

Soft biological tissues are commonly treated as hyperelastic materials for which the consti- 
tutive relations are given in terms of a scalar- valued strain energy-density function *F(F) such 
that 

d¥(F) 

where P is the first Piola-Kirchhoff or the nominal stress tensor. The true or Cauchy stress 
tensor is related to the first Piola-Kirchhoff stress tensor via the "push forward" operation 

a = J l PF T . (5) 

In the absence of body forces the equilibrium equation is 

V • a = 0. (6) 
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The strain energy-density function of an isotropic material can be expressed in terms of the 
three invariants of the Cauchy-Green strain tensor C = F r F, namely 

7i=TrC, / 2 = I(/2-Tr(C 2 )), 7 3 = detC. (7) 

The strain energy-density function of a transversely isotropic (TI) material, whose preferred di- 
rection in the reference configuration is along the unit vector Ni, depends on the two additional 
invariants 

7 4 = Ni-CNi and 7 5 =Ni-C 2 Ni. (8) 

A tissue with two identical families of fibers admits an orthotropic symmetry. If the two fam- 
ilies of fibers are aligned along the unit vectors N\ and N 2 , the strain energy-density function 
describing its behavior depends on three more invariants. Two are analogous to the ones listed 
in Eq. (8), namely 



/ 6 = N 2 -CN 2 and / 7 =N 2 -C 2 N 2 , (9) 



and the last one 



/ 8 = Ni-CN 2 , (10) 

corresponds to the interaction between the two families. 

As was mentioned in the introduction biological tissues are heterogeneous materials. Ac- 
cordingly, we consider next a heterogeneous body made out of N distinct phases whose behav- 
iors are characterized by the strain energy-density functions *p( r ) , r = 1 , . . . , N. In the reference 
configuration each phase occupies a domain 3$<^ r \ The local behavior of the heterogeneous 
material is given in terms of the strain energy-density function 

¥(F,X) = ^(p (r) (X)^W(F), (11) 

where the characteristic functions 

otherwise, 

describe the spatial distributions of the phases within the body. The volume fraction of the 
r-phase in the material is 

C W = I /<pW(X)dV, (13) 



^0 



where V is the referential volume of the body and clearly Y^=i c^ = 1. 

For accurate modeling of this heterogeneous material behavior we need, at one hand to ac- 
count for the interactions between the phases and, at the other, to characterize its macroscopic 
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behavior in terms of relevant information which is elevated from the microscopic level. To ac- 
complish this task we recall the homogenization procedure introduced in the pioneering works 
of Hill (1972) and Ogden (1974). Thus, we apply homogeneous boundary conditions 

x = F X (14) 

on the boundary of the heterogeneous body d£$o, where Fq is a constant matrix with detFo > 0. 
By application of the divergence theorem it can be shown that the average deformation gradient 
satisfies the relation 

F=^|f(X)</V = F . (15) 



The average nominal stress tensor is 



^ J P(X)dV. (16) 



By application of the principle of minimum energy, the effective or macroscopic strain energy- 
density function is (Hill, 1972; Ogden, 1974), 

*P(F) = inf _ { - f ¥(F, X)dV } , (17) 



fgjt(f) V 



$0 



where JT(F) = |f | F = , X e X(X) = FX, X e d^ } is the set of kinematically 
admissible deformation gradients. The corresponding macroscopic constitutive relation is 

The associated macroscopic true stress is determined by the "push forward" operation Eq. (5) 
with respect to F. Generally, application of the variational principle (17) to heterogeneous ma- 
terials can lead to bifurcations corresponding to abrupt changes in the nature of the solution to 
the optimization problem (e.g., Triantafyllidis and Maker, 1985; Rudykh and deBotton, 2011). 
In this work the issue of instabilities is not accounted for and the response of the heterogeneous 
media is considered only along the primary loading curve. 

The solution to the variational statement Eq. (17) is not trivial, and we will make use of the 
NLC variational estimate of deBotton and Shmuel (2010) to obtain close form estimates for the 
effective behavior of the tissue. The NLC variational estimate for *P(F) states that 

❖(F) = ❖o(F) + jV r W [^ r )(F) - ^ r) (F)} , (19) 

r=\ ^ 

where *Pq is an estimate for the effective SEDF of a comparison hyperelastic heterogeneous 

(r) 

material whose phases behaviors are governed by the SEDFs \ and their distributions are 
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characterized by the functions (p^ given in Eq. (12). The term appearing in the second part of 
Eq. (19) is denoted the corrector term, and we note that it depends only on the properties of 
the phases composing the two heterogeneous materials. 

This method was applied by deBotton and Shmuel (2010) and Rudykh and deBotton (201 1) 
to determine the behaviors of transversely isotropic materials with one family of fibers. Thus, it 
is assumed that there are two distinct constituents in the tissue, the fibrous phase and the extra- 
fibrillar matrix. Each phase is incompressible, isotropic, and its behavior can be approximated 
by a hyperelastic model that depends on I\ only. Specifically, we denote by and 
vpO)^) the strain energy-density functions of the fiber and the extra-fibrillar matrix phases, 
respectively. The volume fractions of the two phases are 

c (f) and cW = 1 - c^\ As the 
comparison media the model for neo-Hookean fibrous materials introduced by deBotton et 
al. (2006) is employed (see also deBotton, 2005; Lopez-Pamies and Idiart, 2010). Following 
deBotton and Shmuel (2010), the NLC estimate for the effective strain energy-density function 
is given in terms the minimization problem 

(F) =imn{c«^W (/< m) (I u I 4 , a*)) +c^^ (l[ f) (/i,/ 4 ,»i)) } , (20) 

where 

l[ r) (7iJ 4 ,co 1 )=/4 + ^ + aW(a) 1 )(/ 1 -/4-^), (21) 
yh \ yh J 

I\ and I4 are the invariants of the macroscopic Cauchy-Green strain tensor C = F r F, and 

a^(Q)) = (l-cWa)) 2 and a^ m \co) = (l +c^O)) 2 + c^O) 2 . (22) 

Denote by cb\ the value of C0\ that yields the solution to the minimization problem (20). In some 
cases &\ can be determined analytically otherwise it must be calculated numerically. In gen- 
eral, however, &>\ — &>\ (F). Nonetheless, since the partial derivative of ^ TI ) with respect to C0\ 
identically vanishes at cb\, the expression for the macroscopic stress can be evaluated analyti- 
cally without the need for differentiating CQ\ with respect to F. Consequently, the macroscopic 
nominal stress is 

P (r/) = £ 2c W<)(/^(/ 1 ,/4,«i) 

r=m,f 

+ pF" r , 

(23) 

where (I\ ) is the derivative of *p( r ) (I\ ) with respect to its argument. 



a (r) (Si)F+ (l-a^(ffli)) ( \-L 2 I FN(g)N 



f (r) 
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We model a tissue with two families of fibers as a layered medium made out of transversely 
isotropic fibrous layers such that in the referential configuration the fibers in two adjacent layers 
are alternately directed along the vectors Ni and N2 {e.g., deBotton and Shmuel, 2009). We 
stress that in contrast with artificial composites where, due to the manufacturing process, the 
interaction is between the two layers of fibers, in biological tissues the interaction is between 
the individual fibers belonging to the two families (e.g., Fig. lb). One of the aims of this work is 
to investigate the role of these different interaction hierarchies and the effect of this assumption 
on the predicted overall response. It is assumed that the thickness of the layers is substantially 
smaller than the thickness of the tissue, and that the layers are perfectly bonded to each other. 
In a tissue with two identical families of fibers the unit vectors are symmetrically aligned with 
respect to a symmetry axis (e.g., Xi-axis in Fig. la), and the volume fractions of the two layers 
are identical. The overall symmetry of this two-hierarchies model is orthotropic. 

In a physiological state it is anticipated that the anti-plane shear components are orders of 
magnitude smaller than the in-plane components (e.g., in the coordinate system of Fig. la, F13 
and F23 are substantially smaller than the other components of F). Accordingly, as outlined in 
Appendix A, the overall behavior of the orthotropic tissue is 



(24) 



where I\ 9 14 and 4 are the appropriate orthotropic invariants of C = F r F. Here, (I\,h, (*h) 
and I\ } (I\ ,4 5 G>z) are given in Eq. (21) with 1$ and a>i replacing I4 and (0\ for the second family 
of fibers, respectively. 

We note that the two minimization problems in the expression for ty( TFF ) can be solved 
independently, and denote by (b\ and (&i the values of C0\ and a>i at the minima, respectively. 
Thanks to an argument similar to the one preceding Eq. (23), the following expression for the 
macroscopic nominal stress is derived analytically. 



p{TFF) = £ c (r)^) (^r) { J xJaM ^ [« W )F + (l - a« ((»! )) (l-Z^Wx®^ 

r=m.f L V / 

+ £ cW^ r) (i< r) (7i,^dfe)) [aW(dfe)F+(l-aW((%))ri-7 6 ^W 2 ®S 2 



r=m,f 

+ pF~ T , 



(25) 
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where, as before, is the derivative of with respect to its argument. We em- 

phasize that this is a close-form explicit expression for the stress and the values for cb\ and 
(i>2 can be substituted once the minimization problem Eq. (24) is solved, either analytically or 
numerically. 

Before we proceed we emphasize that the expressions for the macroscopic SEDF *y( TFF ) 
and the macroscopic stress p( TFF ) involve, as they should, the histology of the tissue in terms 
of the direction and volume fraction of the collagen fibers together with the expressions for the 
behaviors of the two constituents in terms of their SEDFs *p( m ) and x ¥^\ 

We consider next the FE model to which the aforementioned analytical expressions are com- 
pared. Drawings of the elementary unit cell and a representative section of the periodic me- 
dia containing four elementary cells together with the chosen coordinate system are shown in 
Fig. la and b, respectively. Arbitrarily we define a material coordinate system such that in the 
reference state the direction between the two families is aligned with the Xi-axis. With respect 
to an artery, axes X\ and X2 may be associated with the circumferential and the axial directions, 
respectively. 

The referential angle between the fiber's direction and the X\ -axis is denoted 0. Accordingly, 
the ratio between the lengths of the planar edges of the cell is 

-=tan0, (26) 
a 

where a and b are the lengths of the faces along the X\ and the X2 axes, respectively. Due to 
the periodicity, the total length of the fiber segments in the unit cell is 2Va 2 + b 2 , and the total 
volume of the fiber in the cell is 

V f = *d 2 Va 2 + b 2 , (27) 
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where d is the diameter of the fibers. Once the volume fraction of the fiber is set, the 
condition for the distance between the planes of the two families h (half the hight of the unit 
cell) is 

' A " ■-■ (28) 



d ) 4cW sin© \a 
The geometrical requirements h/d > 1 and a sin©/ d > 1 imply that 

1 < \ < ^77Y- (29) 
d 4c(/) 

We note that in the (Xi,^) -plane the spacing between the fibers is a sin©, and the spacing 
between the fibers in the X?> direction is h. A requirement that the spacing will be identical 
leads to 

" ^ (30) 



(32) 



d V 4c(/) 

The response of the heterogeneous material is obtained by applying periodic boundary con- 
ditions on a single unit cell {e.g., Fig. la). In the undeformed configuration the unit cell occu- 
pies the domain 

(KXi^a, (KX 2 0, 0^X 3 ^2h. (31) 

The expressions for these conditions are written in terms of the components of the macroscopic 
deformation gradient tensor F as follows: 

(1) The displacements of points on the right (X\ — a) and the left (X\ = 0) faces are related 
via 

u\ = u^+ (F\i — l^ja 

u\ — U2+F21CL 

and traction free conditions were assumed for the shear components in the X3 direction. 

(2) The displacements of points on the front (X2 = 0) and the rear (X2 = b) faces are related 
via 

u F , = u f E +F n b 

/ x ' (33) 
u F 2 =uf +[F 22 -l)b 

and traction free conditions were assumed for the shear components in the X3 direction. 

(3) On the top (X3 = 2h) and bottom (X3 = 0) faces = u F = 0, and traction free condi- 
tions were assumed for the shear components in the X\ and X 2 directions. 

The overall behavior of the heterogeneous media was then determined from the relations be- 
tween the average of the stress field (e.g., Eq. (16) with Eq. (5)), and the corresponding applied 
macroscopic deformation F. 
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For inclined loadings we distinguish, for convenience, between the principal coordinate sys- 
tem of the right Cauchy-Green deformation tensor and the material coordinate system. Under 
plane-strain loading condition, thanks to the assumed incompressibility of the tissue, in the 
principal coordinate system the most general macroscopic deformation gradient is 

( X \ 



F' 



(34) 



/ 



X~ l 
\ 1 

where X is the principal stretch ratio. This is related to the average deformation gradient in the 
material coordinate system via the relation 



R r F'R, 



(35) 



where 



/ cos/3 -sin/3 \ 
R= sin/3 cos/3 , (36) 
V 1/ 

and j8 is the referential angle between the direction of the principal stretch ratio and the axis 
Xi. 

3. Applications 

In this section we make use of both the variational estimate and the periodic FE simulation 
technique to determine the overall behavior of a specific tissue with two families of fibers. 
To this end we first make an appropriate choice for the SEDFs of the collagen fibers and the 
extra-fibrillar matrix. 

While most of the biochemical processes in the tissue occur within the extra-fibrillar matrix, 
from a mechanical viewpoint the fiber phase is more interesting. In the referential state of the 
tissue the fibers are crimped and straighten when the tissue is subjected to tension in the fiber 
direction (e.g., Freed and Doehring, 2005; Holzapfel, 2008). Thus, initially the contribution of 
the fiber is quite small but when straighten it becomes the primary load carrying component in 
the tissue. Since not all the fibers straighten simultaneously, there is a recruitment stage during 
which more and more fibers straighten and gradually increase the overall stiffness of the tissue. 
A hyperelastic model that can capture this phenomenon is the Gent (1996) model 

/i-3> 



¥ G (F;jMm) = ~tiJ m ^ (l 



(37) 



Here \i is a shear like modulus and J m is a dimensionless "locking" parameter. Qualitatively, 
the behavior exhibited by this model can be characterized as follows. In the small strains regime 
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when I\ <C J m + 3, the linear term of Taylor series expansion of Eq. (37), namely ^jih, is the 
dominating term. In accordance with the recruitment processes, as I\ increases the material 
stiffens. Finally, in the limit (h — 3)—>J m there is a dramatic stiffening and the material locks 
up as becomes unbounded. For example, in the pure shear loading condition of Eq. (34) 



According to the findings of Gundiah et al. (2007) a hyperelastic model that can capture the 
behavior of the extra-fibrillar matrix with the elastin fibers is the incompressible neo-Hookean 
model 



This model is characterized by a single shear-like modulus ji. In the limit J m —> °o the Gent 
model (37) reduces to the neo-Hookean model. 

In the course of this work models of tissues with neo-Hookean matrices were analyzed by 
both the NLC variational estimate and the finite element method, and compared with corre- 
sponding predictions of phenomenological models. It was found that as long as no dispersion 
of the fibers is assumed the usage of a neo-Hookean model for the matrix phase leads to a fine 
agreement with the phenomenological predictions. However, as soon as the distributions of the 
fibers about their primary directions are accounted for, the neo-Hookean model fails to provide 
satisfactory agreement with the experimentally fitted findings. The reason is the inability of 
the neo-Hookean model to capture the stiffening of the dispersed fibers during extensions that 
are not along the primary directions. Consequently, in the current study that is aimed towards 
demonstration of the ability of the micromechanics based model to capture the essential as- 
pects of the behavior of heterogeneous soft materials, we assume no dispersion of the fiber and 
a Gent model for the behavior of the extra-fibrillar matrix. The neo-Hookean model will be 
used in subsequent studies where dispersion of the fibers will be accounted for. 

Aside from the ability of the chosen Gent model to capture the behaviors of the two constitut- 
ing phases, a computational advantage stems from the fact that with this choice the optimization 
problem (24) can be solved analytically. Thus, the associated Euler-Lagrange equations admit 
the form of cubic polynomials in C0\ and o>2, from which explicit close-form expressions for 
(b\ and (&i are obtained (deBotton and Shmuel, 2010). We emphasize, however, that within the 
framework of the proposed variational estimate (24), the models for the tissue constituents can 
be easily modified or replaced with other, possibly more sophisticated, models. 

In total there are 4 material parameters corresponding to the shear moduli and the locking 
parameters of the extra-fibrillar matrix and the fibers. In addition, there are two histological 
parameters, namely the volume fraction of the fibers and the opening angle between the two 



the critical stretch ratio at which locking occurs is X c = ^ (V^~+4+ y/Jm) • 




(38) 
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FIGURE 2. Stored elastic energy-density as a function of the stretch ratio according to 
NLC estimate (continuous curves), the periodic FE model (shaded marks), and the model of 
Holzapfel et al. (2005) (dashed curves). Figures (a) and (b) correspond to extensions along the 
in-plane principal orthotopic axes X\ and X2, respectively. 



families of fibers. The parameters were chosen by fitting the overall prediction of the microme- 
chanics based models to the predictions of a phenomenological model whose own parameters 
were fitted to experimental data. Variety of results are available in the literature among which 
we recall the works of Holzapfel (2006); Hariton et al. (2007a); Kiousis et al. (2009); Stalhand 
(2009) and Mortier et al. (2010). Here we calibrated the micromechanics model with the aid 
of the strain energy-density function given in Holzapfel et al. (2005) that accounts for the dis- 
persion of the fibers. This model involves five materials parameters. Mean values for these 
parameters are also given in Holzapfel et al. (2005) for the three layers composing the wall 
of healthy human coronary arteries. From a mechanical viewpoint the relevant layers are the 
Adventitia and the Media (Holzapfel, 2008). At this stage of our study we wish to treat the 
tissue as a single layer and hence we used the average value of the parameters for the two 
dominant layers. Thus, the numerical values used for the five parameters appearing in Eq. (1) 
of Holzapfel et al. (2005) are \i = 4.42 kPa, k x = 30.1kPa, k 2 = 46.6, p = 0.4 and half the 
opening angle between the two families of fibers = 43°. Note that these values are also in 
fair agreement with corresponding results that were reported in Stalhand (2009). Additionally, 
Holzapfel et al. (2000); Balzani et al. (2006); Holzapfel (2006) and Hariton et al. (2007a) re- 
ported quite similar estimates for the angle between the two fiber families. The estimate for the 
collagen volume fraction was deduced from the experimental findings of de Figueiredo Borges 
et al. (2008) and the fitting process of deBotton and Shmuel (2009). 
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FIGURE 3 . The deformed configuration of the unit cell during extensions alon (a) the X\ and 
(b) the X2 directions. The stored energy-density at several slices is depicted by the contour plots. 

The numerical values that were chosen for the modeled tissue are summarized in table 1. 
These were determined by comparing the predictions for the stored energy-density according 
to the micromechanics based model Eq. (24) and the FE periodic model with the model of 
Holzapfel et al. (2005) during extensions along the principal orthotropic axes. This comparison 
is depicted in Fig. 2. The continuous curves depict the predictions of the NLC estimate, the 
dark circle and diamond marks those of the periodic FE model, and the dashed curves to the 
model of Holzapfel et al. (2005). Figs. 2(a) and (b) correspond to extensions along the in- 
plane principal orthotropic axes X\ and X2, respectively. We note that this choice of parameters 
leads to a fine agreement between the three different models and emphasize that throughout 
this work the parameters listed in Table 1 are used for both the NLC variational estimate and 
the FE simulation. We recall that for the periodic FE model an additional parameter, namely 
the relative spacing between the fibers, needs to be set. For the chosen histological parameters 
inequality (29) is 1 < h/d < 2.6, and according to Eq. (30) the choice hjd — 1.62 will lead 
to identical spacing. We note that a choice of hjd < 1.62 will increase the cross-families 
interaction between the fibers in the two families, whereas a reversed choice will increase the 
internal interaction between the fibers in the two families. In this work the choice hjd — 1.2 
was made in order to maximize the interaction between the fibers in the two families. This was 
done in order to critically examine the assumption embedded in the analytical model that the 
out of plane interaction is between layers of fibers and not between the fibers themselves. 



Table 1. Material parameters 





Jm) 




Af) 

J m 


c (f) 





3.5 kPa 


0.098 


280 kPa 


0.45 


0.30 


43° 
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Min: 194.238 C^O Min: 154 ' 621 

FIGURE 4. Deformed configuration of the unit cell during incline extension at /3 = 21.5°. 
Case-(/) is shown in the left figure and case-(n) in the right one. The stored energy-density at 
several slices is depicted by the contour plots. 



The deformations of the unit cell and the distributions of the stored strain energy during 
the two loading conditions are demonstrated in Fig. 3 for X — 1.29. In both cases the fibers 
are stretched and the elastic energy stored in the fibers is an order of magnitude higher than 
the one in the EFM. This is due to the fact that the stiffer fibers carry most of the mechanical 
load at this stage of the deformation. It is interesting to note that although the stretching is 
along the principal orthotropic axes, the cell undergoes local shear and hence in the deformed 
configuration its boundary is not straight. This local shear is associated with the imposed 
periodic boundary conditions. 

To examine the ability of the NLC estimate to capture the behavior predicted by the periodic 
FE simulation, we also compare their predictions under incline loading conditions. Two con- 
figurations are considered corresponding to tensions at angles j8 = 21.5° and j8 = 43° relative 
to the principal orthotropic Xi-axis. These loading states may correspond to combinations of 
axial and hoop extensions together with a twist of the artery. Particularly, with j8 = 43° the ten- 
sion is along the first fiber families (the bottom family in Fig. 1), while the second family is in a 
state of compression in the fiber direction. Under tension at an angle j8 = 21 .5°, the first family 
is under a combination of shear and tension while the second family is under compression and 
shear. 

It is widely accepted that the collagen fibers do not contribute to the stiffness of the tissue 
when compressed in the fiber direction (e.g., Holzapfel, 2008). Accordingly, during the incline 
tests we examined two different cases. In case-(/) the behavior of the compressed fiber family 
was set to be identical to that of the extra-fibrillar matrix. That is, under compression the fiber 
behaves according to the Gent model with ji^ and Essentially, this resulted in models 
with a single load-carrying family of fibers, the one under tension. In case-(//) the fiber behavior 
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FIGURE 5. Deformed configuration of the unit cell during incline extension at j8 = 43°. 
Case-(/) is shown in the left figure and case-(//) in the right one. The stored energy-density at 
several slices is depicted by the contour plots. 

was not modified and the contribution of the second family to compression was accounted for 
(i.e., the fiber behaves according to the Gent model with ji^ and jffl)- 

Deformed configurations of the periodic cell together with contour plots representing the 
stored energy are shown in Figs. 4 and 5 for the inclined extensions at j8 = 21.5° and j8 = 43°, 
respectively. The severe macroscopic shear to which the cell is subjected can be appreciated. 
In addition, the local shear resulting from the periodicity conditions is associated with the 
curvature of the cell faces. In case-(ii) the upper family of fibers which is subjected to com- 
pression can be easily distinguished due to the higher level of stored energy in comparison with 
the surrounding matrix. Contrarily, in case-(/) the amount of the stored energy-density in the 
compressed fibers is identical to the one in the surrounding EFM. 

In Fig. 6 the estimates for the stored energy-density according to the two micromechanics 
based models are shown. Figs. 6(a) and (b) correspond to extensions at angles j8 = 21.5° and 
j8 = 43°, respectively. The NLC estimate is represented by the continuous curves, the periodic 
FE model by the dark marks, and we also depict the corresponding results according to the 
model of Holzapfel et al. (2005) (dashed curves). Case-(/), with no contribution of the fibers 
in compression, is identified by blue square marks and case-(//), where their contribution is 
accounted for, by the green diamond marks. 

Our first observation concerns the excellent agreement between the two micromechanics 
based models. Thus, in all four cases the shaded marks that correspond to the FE simulations 
are lying on top of the continuous curves derived by the NLC estimate. This suggest that 
the close form NLC estimate accurately captures the response of the heterogeneous material. 
Contrarily, while the general trend of the phenomenological model is reminiscent of the mi- 
cromechanics models, the precise details are different. This is because the different structure of 



ANALYSIS OF THE MICROMECHANICS OF SOFT FIBROUS CONNECTIVE TISSUES 



18 



-a 
i 

>^ 

a 
-a 
o 

+-> 

GO 



(a) 




a 

i 

to/) 

id 

o 
55 



(b) 



1 T/ 3 

kJ/m 






1 1 
/ 1 










1 1 

i 

i [ 
if 

Ui 










fi i 
? 1 1 

/' i r 
J' i J 








V < 

/ ' J> 


1 1 r 
i J 
if 











1.05 



1.1 



1.15 



FIGURE 6. The stored elastic energy-density according to the NLC estimate (continuous 
curves), the periodic FE model (shaded marks), and the model of Holzapfel et al. (2005) (dashed 
curves). Figures (a) and (b) correspond to tensions in the j3 = 21.5° and the j3 = 43° directions, 
respectively. Cases (/) and (ii) are signified by the blue square and the the green diamond marks, 

respectively. 



1.2 



the latter, and its inability to capture the local interactions between the collagen fibers and their 
surrounding. Particularly, the excessively stiff response predicted by the phenomenological 
model is due to the inherently embedded assumption that the fibers and the tissue are subjected 
to the same macroscopic deformation. 

An interesting observation concerns the influence of the behavior of the fibers under com- 
pression. The difference between cases (/) and (ii) is represented by the difference between 
the curves marked by blue squares and those marked by green diamonds, respectively. We 
also computed the stored energy-density according to the phenomenological model for these 
two cases. Ones according to case-(/) without the contribution of the compressed family (blue 
dashed curve with hollow square marks), and once according to case-(//) with the contribution 
of this family (green dashed curve with hollow diamond marks). 

As anticipated, the curves corresponding to case-(/) are laying beneath the ones for case-(//). 
However, while initially the difference between the two pairs of curves increases, eventually 
the two are getting closer. The reason is the rotation of the fibers towards the tension direction 
during the deformation. Because of this, at some value of X the compressive strains in the 
compressed family start to decrease till at large enough stretch ratio they turn into tension 
(deBotton and Shmuel, 2009). The fact that the level of the stored energy in the compressed 
family is severely lower than the one in the tensed family can be seen in Figs. 4(ii) and 5(h)- It 
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FIGURE 7 . Evolution of the angle between the fiber directions and the principal orthotropic 
axis versus the tensile stretch ratio. The curves correspond to the analytical NLC estimate and 
the marks to the periodic FE model. 

is important to note that according to the micromechanics based models the difference between 
the two cases is larger than the difference predicted by the phenomenological model. 

The rotation of the fibers, in terms of the evolution of the angle between the fiber and the 
principal orthotropic direction X\ versus the principal stretch ratio is shown in Fig. 7 for the 
four different loadings shown in Figs. 3, 4 and 5. The continuous curves correspond to the 
results of the NLC estimate and the marks to the FE simulations. Following the results of 
Appendix A, the analytical expression for the angle is simply 

9 = arctan (^^ +/J . (39) 

There is an excellent agreement between the two micromechanics based models. 

To assess the ability of the NLC estimate to capture the behavior of the heterogeneous tissue 
we also compare the overall stresses determined via Eq. (25) with the corresponding stresses 
predicted by the finite element simulation according to Eq. (16). In Fig. 8 the predicted compo- 
nents of the deviatoric stress according to the different models are shown during tensions along 
the principal orthotropic axes of the tissue. Figs. 8a and 8b correspond to tensions along the 
in-plane orthotropic axes X\ andX2, respectively. The corresponding deformed configurations 
are shown in Figs. 3a and 3b, respectively. Throughout, the curves correspond to the predic- 
tions of the NLC estimate, and the marks to the FE simulations. The blue curves and square 
marks correspond to the stress in the tensile direction (e.g., in the coordinate system of Fig. 1, 
of[ for Fig. 8a and <7^ for Fig. 8b). The compressive stresses in the perpendicular direction 
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(a) A (b) A 

FIGURE 8. Deviatoric stresses in the tension direction (blue squares), in the transverse di- 
rection (green circles), and in-plane shear (red diamonds) versus the tensile stretch ratio. The 
continuous curves correspond to the NLC estimate and the marks to the periodic FE model. 
Figures (a) and (b) depict tensions along the X\ and X2 axes, respectively, respectively. 

are represented by the green curves and the circular marks. The in-plane shear stresses are 
represented by the red curves and the diamond marks. 

The overall response of the tissue, which is dictated by the behavior of the collagen fiber, can 
be traced in these figures. Initially, the contribution of the fiber is negligible and the slope of 
the stress-strain curve is very low. At a stretch ratio of about 1.2 the contribution of the fibers 
becomes evident as the intensity of the stresses become larger than the shear moduli of the 
EFM. Subsequently, as the stretch ratio approaches 1.3, the stiff fibers dominate the response of 
the tissue and the slopes of the stress-strain curves increase dramatically. In both figures, since 
the stretches are along the principal orthotropic axes, the macroscopic shear stresses vanish. 
In the context of the present work, we emphasize the fine agreement between the analytical 
prediction according to the NLC estimate and the predictions according to the FE simulations. 

The stresses developing due to incline extensions are shown in Fig. 9. Figs. 9a and 9b 
correspond to tensions at j8 = 21.5° and j8 =43°, respectively. The corresponding deformed 
configurations are shown in Figs. 4 and 5, respectively. Throughout, the curves correspond to 
the predictions of the NLC estimate, and the marks to the FE simulations. Case-(/) is repre- 
sented by the continuous curves and the square, circle and diamond marks, and Case-(//) by 
the dashed curves and the upright triangle, inverted triangle and cross marks. The blue curves, 
square and upright triangle marks correspond to the stress in the tensile direction. The per- 
pendicular compressive stresses are represented by the green curves, circle and the inverted 
triangle marks. The in-plane shear stresses are represented by the red curves, the diamond and 
the cross marks. 
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FIGURE 9 . Deviatoric stresses in the tension direction (blue squares and upright triangles), in 
the transverse direction (green circles and inverted triangles), and in-plane shear (red diamonds 
and crosses) versus the tensile stretch ratio. The curves correspond to the NLC estimate and the 
marks to the periodic FE model. Case-(z) corresponds to the continuous curves and the square, 
circle and diamond marks. Case-(z7) corresponds to the dashed curves and the upright triangle, 
inverted triangle and cross marks. 



Once again, we observe that throughout the marks are laying on top of the curves. This 
indicates that there is an excellent agreement between the analytical NLC estimate and the cor- 
responding FE simulations. The role of the fibers can be appreciated with regard to the overall 
response of the tissue. In the case of tension at j8 = 21.5°, since the direction of the tension 
is closer to the direction of one of the families (the lower family in Fig. 1), the contribution of 
the fiber becomes evident at a faster rate and already at a stretch ratio of A = 1.15 the tissue 
stiffens. Naturally, during a tension along one of the families (j8 = 43°) this phenomenon is 
further enhanced and the tissue stiffens already at X — 1.1. An interesting observation con- 
cerns the shear stresses. When the extension is at j8 = 21.5° large shear stresses develop in the 
tissue since the principal stretch is shifted from the principal orthotropic axes. At j8 = 43° the 
shear stress are much smaller since the load is carried by a single family and the rotation of 
the second family is negligible (see also Fig. 7). In particular, since in case-(/) the extension is 
essentially along the preferred direction of the tissue, no shear stresses develop. 



4. Concluding remarks 

Highly accurate and efficient methods for characterizing the mechanical response of biolog- 
ical tissues are required in modern analyses and modeling. Micromechanics based approaches 
that enable to account for the behaviors of the individual constituents together with the histol- 
ogy of the tissue can provide realistic and versatile methods for characterizing the behaviors of 
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broad classes of tissues. For instance, if estimates for the behaviors of the collagen fibers and 
the EFM are available, by following this approach different parts of the arterial wall, different 
sections of the artery, and even different arteries can be modeled by modifying appropriate 
morphological parameters such as the volumetric content and the spatial arrangement of the 
fibers. Moreover, it is possible that many of the histological parameters will be measured by 
non-invasive techniques, and hence the need for large sets of experiments required for tuning 
up phenomenological models can be avoided. 

In this work we made use of micromechanics based estimates for soft connective tissues with 
one and two families of fibers. Both analytical and numerical models were examined and their 
predictions for the overall response of a healthy human coronary artery were compared. The 
analytical model, which is based on the variational estimate of deBotton and Shmuel (2010), 
resulted in a close form expression which is no more complicated than corresponding ad-hoc 
phenomenological models. The numerical model is based on the homogenization of a periodic 
unit cell subjected to periodic boundary conditions. The constitutive behaviors for the collagen 
fibers and the EFM were calibrated by comparison with a phenomenological model whose pa- 
rameters were determined for the modeled tissue during circumferential and axial tensile tests. 
Under inclined loadings, when one of the families is tensed and the other is compressed, we 
examined the behavior of the tissue when the compressed family either contributes to the over- 
all response or not. We find that under inclined loadings there is marked difference between 
the predictions of the phenomenological and the micromechanical estimates, where the former 
predicts a substantially stiffer response of the arterial wall. Contrarily, for all loading condi- 
tions there is a fine agreement between the predictions for the overall stresses and evolution 
according to the two micromechanics based models. This implies that the proposed close form 
analytical estimate adequately captures the response of connective tissues and the interactions 
between the EFM and the two families of collagen fibers. 

Appendix A. The effective behavior of an incompressible laminated medium 

We follow the steps followed by deBotton (2005) and determine the effective strain energy- 
density function *p( L ) (F) of a layered medium made out of alternating incompressible layers 
whose behaviors are characterized by the strain energy-density functions (F) , / = 1 , 2. We 
assume that the thickness of the individual layers is at least an order of magnitude smaller 
than the thickness of the medium under consideration which, in turn, is smaller than its other 
dimensions. This scale separation assumption implies that under homogeneous boundary con- 
dition of Eq. (14) the local deformation gradients in the layers F^ and F( 2 ) are constants. 
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Accordingly, from Eq. (15) 

F= £ c (/) F (/) = F , (40) 
l=\2 

where are the volume fractions of the layers composing the medium such that Y*i=i 2 c ^ — 
1. The continuity of the displacements across the interfaces between two adjacent layers are 

(fW-F (2) )M l = 0, (41) 

where is any unit vector in the layers plane such that Ml • = with being the 
unit vector normal to the layers plane. Thus, only the components of the normal projections 
of the local deformation gradients F^N^ may differ. Specifically, in the coordinate system of 
Fig. la, only the components F13, F23 and F33 may differ. If we assume that the anti-plane shear 
components are vanishingly small (i.e., (FNl) • Ml = for any vector M^, or in the chosen 
coordinate system F13 = F23 = 0), then in each layer (F^N^) • Nl = (FNl) • Nl since due to 
the incompressibility constraint these depend on the continuous in-plane components (in the 
chosen coordinate system — F33). In summary, Eqs. (40), (41) and the incompressibility 
assumption lead to the conclusion that under vanishingly small anti-plane shear loads 

FW =F, /= 1,2. (42) 

Consequently, the solution of the optimization problem Eq. (17) is trivial, and the effective 
energy-density function of the layered medium is nothing but the volume average of the con- 
stituents energy-density functions, namely 

xpW(F)= £ S¥ l \F). (43) 

Z=l,2 

We finally note that the above developments are independent of the precise constitutive behav- 
iors of the phases composing the laminate. Moreover, the above arguments can be trivially 
extended to laminated media with more than two distinct phases (i.e., I > 2). This last obser- 
vation may become useful if a distribution of the fibers is accounted for. 

References 

Agoras, M., Lopez-Pamies, O. and Castaneda, P. P. (2009). A general hyperelastic model for incom- 
pressible fiber-reinforced elastomers, /. Mech. Phys. Solids 57: 268-286. 

Alastrue, V., Martinez, M. A., Doblare, M. and Menzel, A. (2009). Anisotropic micro-sphere-based 
finite elasticity applied to blood vessel modelling, /. Mech. Phys. Solids 57: 178-203. 

Ambrosi, D., Ateshian, G. A., Arruda, E. M., Cowin, S. C, Dumais, J., Goriely, A., Holzapfel, G. A., 
Humphrey, J. D., Kemkemer, R., Kuhl, E., Olberding, J. E., Taber, L. A. and Garikipati, K. (2011). 
Perspectives on biological growth and remodeling, /. Mech. Phys. Solids 59: 863-883. 



ANALYSIS OF THE MICROMECHANICS OF SOFT FIBROUS CONNECTIVE TISSUES 



24 



Avery, N. C. and Bailey, A. J. (2008). Restraining cross-links responsible for the mechanical properties 
of collagen fibers: natural and artificial, in P. Fratzl (ed.), Collagen Structure and Mechanics, Springer 
Science Business Media, Inc., New York, chapter 4. 

Balzani, D., Neff, P., Schroder, J. and Holzapfel, G. A. (2006). A polyconvex framework for soft 
biological tissues. Adjustment to experimental data, Int. J. Solids Structures 43: 6052 - 6070. 

Brun, M., Lopez-Pamies, O. and Ponte Castaneda, P. (2007). Homogenization estimates for fiber- 
reinforced elastomers with periodic microstructures, Int. J. Solids Struct. 44: 5953-5979. 

Chen, H., Liu, Y., Zhao, X., Lanir, Y. and Kassab, G. S. (2011). A micromechanics finite-strain consti- 
tutive model of fibrous tissue, /. Mech. Phys. Solids 59: 1823-1837. 

de Figueiredo Borges, L., Jaldin, R. G., Dias, R. R., Stolf, N. A. G., Michel, J. B. and Gutierrez, P. S. 
(2008). Collagen is reduced and disrupted in human aneurysms and dissections of ascending aorta, 
Human Pathology 39: 437-443. 

deBotton, G. (2005). Transversely isotropic sequentially laminated composites in finite elasticity, /. 
Mech. Phys. Solids 53: 1334-1361. 

deBotton, G. and Hariton, I. (2006). Out-of-plane shear deformation of a neo-Hookean fiber composite, 
Physics Letters A 354: 156-160. 

deBotton, G. and Shmuel, G. (2009). Mechanics of composites with two families of finitely extensible 
fibers undergoing large deformations, /. Mech. Phys. Solids 57: 1165-1181. 

deBotton, G. and Shmuel, G. (2010). A new variational estimate for the effective response of hypere- 
lastic composites, /. Mech. Phys. Solids 58: 466-483. 

deBotton, G., Hariton, I. and Socolsky, E. A. (2006). Neo-Hookean fiber-reinforced composites in finite 
elasticity, /. Mech. Phys. Solids 54: 533-559. 

Eberth, J. E, Popovic, N., Gresham, V. C, Wilson, E. and Humphrey, J. D. (2010). Time course of 
carotid artery growth and remodeling in response to altered pulsatility, Am. J. Physiol. Heart Circ. 
Physiol. 299: H1875-H1883. 

Fratzl, P. (ed.) (2008). Collagen Structure and Mechanics, Springer Science Business Media, Inc., New 
York. 

Freed, A. D. and Doehring, T. C. (2005). Elastic model for crimped collagen fibrils, /. Biomech. Eng., 

Trans. ASME 127: 587-593. 
Fung, Y. C. (1967). Elasticity of soft tissues in simple elongation, Am. J. Physiol. 28: 1532-1544. 
Gent, A. N. (1996). A new constitutive relation for rubber, Rubber Chemistry and Technology 69: 59-61 . 
Gundiah, N., Ratcliffe, M. B. and Pruitt, L. A. (2007). Determination of strain energy function for 

arterial elastin: Experiments using histology and mechanical tests, Journal of Biomechanics 40: 586- 

594. 

Hariton, I., deBotton, G., Gasser, T. C. and Holzapfel, G. A. (2007a). Stress-driven collagen fiber 
remodeling in arterial walls, Biomechanics and Modeling in Mechanobiology 6: 163-175. 

Hariton, I., deBotton, G., Gasser, T. C. and Holzapfel, G. A. (2007b). Stress modulated collagen fibers 
remodeling in a human carotid bifurcation, /. Theor. Biol. 248: 460-470. 

Hill, R. (1972). On constitutive macro- variables for heterogeneous solids at finite strain, Proc. R. Soc. 
Lond. A 326: 131-147. 

Holzapfel, G. A. (2006). Determination of material models for arterial walls from uniaxial extension 

tests and histological structure, /. Theor. Biol. 238: 290-302. 
Holzapfel, G. A. (2008). Collagen in arterial walls: Biomechanical aspects, in P. Fratzl (ed.), Collagen 

Structure and Mechanics, Springer Science Business Media, Inc., New York. 



ANALYSIS OF THE MICROMECHANICS OF SOFT FIBROUS CONNECTIVE TISSUES 



25 



Holzapfel, G. A. and Ogden, R. W. (2010). Constitutive modelling of arteries, Proc. R. Soc. Lond. A 
466: 1551-1597. 

Holzapfel, G. A., Gasser, T. C. and Ogden, R. W. (2000). A new constitutive framework for arterial wall 

mechanics and a comparative study of material models, /. of Elasticity 61: 1-48. 
Holzapfel, G. A., Sommer, G., Gasser, C. T. and Regitnig, R (2005). Determination of layer- specific 

mechanical properties of human coronary arteries with nonatherosclerotic intimal thickening and 

related constitutive modeling, Am. J. Physiol. Heart Circ. Physiol. 289: H2048-H2058. 
Holzapfel, G. A., Stadler, M. and A., S.-B. C. (2002). A layer- specific three-dimensional model for the 

simulation of balloon angioplasty using magnetic resonance imaging and mechanical testing, Annals 

of Biomedical Engineering 30: 753-767 '. 
Humphrey, J. (2003). Review paper: Continuum biomechanics of soft biological tissues, Proceedings of 

the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences 459: 3-46. 
Humphrey, J. D. (2002). Cardiovascular Solid Mechanics: cells, tissues, and organs, Springer- Verlag, 

New York. 

Humphrey, J. D. (2008). Vascular adaptation and mechanical homeostasis at tissue, cellular, and sub- 
cellular levels, Cell Biochem. Biophys. 50: 53-78. 

Humphrey, J. D. (2009). Vascular mechanics, mechanobiology, and remodeling, Journal of Mechanics 
in Medicine and Biology 9: 243-257. 

Humphrey, J. D. and Yin, F. C. (1987). A new constitutive formulation for characterizing the mechanical 
behavior of soft tissues, Biophysical Journal 52: 563 — 570. 

Kiousis, D. E., Rubinigg, S. R, Auer, M. and Holzapfel, G. A. (2009). A methodology to analyze 
changes in lipid core and calcification onto fibrous cap vulnerability: The human atherosclerotic 
carotid bifurcation as an illustratory example, /. Biomech. Eng., Trans. ASME 131(12): 121002. 

Lanir, Y. (1983). Constitutive equations for fibrous connective tissues, /. Biomech. Eng., Trans. ASME 
16: 1-12. 

Lopez-Pamies, O. and Idiart, M. (2010). Fiber-reinforced hyperelastic solids: a realizable homogeniza- 
tion constitutive theory, /. of Engineering Mathematics 68: 57-83. 

Menzel, A. and Waffenschmidt, T. (2009). A microsphere-based remodelling formulation for 
anisotropic biological tissues, Phil. Trans. R. Soc. A: Mathematical, Physical and Engineering Sci- 
ences 367: 3499-3523. 

Mortier, R, Holzapfel, G., De Beule, M., Van Loo, D., Taeymans, Y, Segers, R, Verdonck, R and 
Verhegghe, B. (2010). A novel simulation strategy for stent insertion and deployment in curved 
coronary bifurcations: Comparison of three drug-eluting stents, Annals of Biomedical Engineering 
38: 88-99. 

Nerem, R. M. and Seliktar, D. (2001). Vascular tissue engineering, Annual Review of Biomedical Engi- 
neering 3: 225-243. 

Ogden, R. W. (1974). On the overall moduli of non-linear elastic composite materials, /. Mech. Phys. 
Solids 22: 541-553. 

Ponte Castaneda, P. and Tiberio, E. (2000). A second-order homogenization method in finite elasticity 
and applications to black-filled elastomers, /. Mech. Phys. Solids 48: 1389-1411. 

Rudykh, S. and deBotton, G. (2011). Instabilities of hyperelastic fiber composites: micromechanical 
versus numerical analyses, /. of Elasticity 106: 123-147. 

Shmuel, G. and deBotton, G. (2010). Out-of-plane shear of fiber composites at moderate stretch levels, 
/. of Engineering Mathematics 68: 85-97. 



ANALYSIS OF THE MICROMECHANICS OF SOFT FIBROUS CONNECTIVE TISSUES 



26 



Stalhand, J. (2009). Determination of human arterial wall parameters from clinical data, Biomechanics 
and Modeling in Mechanobiology 8: 141-148. 

Taber, L. A. and Humphrey, J. D. (2001). Stress-modulated growth, residual stress, and vascular hetero- 
geneity, /. Biomech. Eng., Trans. ASME 123: 528-535. 

Triantafyllidis, N. and Maker, B. N. (1985). On the comparison between microscopic and macro- 
scopic instability mechanisms in a class of fiber-reinforced composites, /. Appl. Mech., Trans. ASME 
52: 794-800. 



